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Abstract 

We consider a spin- 1/2 charged particle in the plane under the influence 
of two idealized Aharonov-Bohm fluxes. We show that the Pauli operator 
as a differential operator is defined by appropriate boundary conditions at 
the two vortices. Further we explicitly construct a basis in the deficiency 
subspaces of the symmetric operator obtained by restricting the domain 
to functions with supports separated from the vortices. This construction 
makes it possible to apply the Krein's formula to the Pauli operator. 
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I. Introduction 



The goal of the present paper is to provide a more detailed analysis of the Pauli 
operator describing a spin- 1/2 charged particle under the influence of two Aharonov- 
Bohm (AB) fluxesP We consider the idealized setup when the magnetic fluxes are 
concentrated along two parallel lines so that the problem effectively reduces to a two- 
dimensional quantum system living in a perpendicular plane. In what follows we call 
the intersection points of the fluxes with the plane vortices. 

To define the Pauli Hamiltonian with singular fluxes we use the Aharonov-Casher 
decompositionP It makes it possible to introduce the two diagonal components of the 
Pauli operator corresponding to spin up and down as the unique selfadjoint operators 
associated to appropriate quadratic forms. Since the magnetic field vanishes outside 
of the vortices the two components of the Pauli Hamiltonian as well as the spinless 
AB Hamiltonian are selfadjoint extensions of the same symmetric operator. In the 
case of one AB flux all the selfadjoint extensions are known to be defined by appro- 
priate boundary conditions at the vortexP^I Thus our first goal was to distinguish the 
boundary conditions defining the two components of the Pauli Hamiltonian. 

The second goal was to construct a basis in the deficiency subspaces in the two- 
vortex case. In this well the two diagonal components of the Pauli Hamiltonian 
and the spinless Hamiltonian are selfadjoint extensions of a common symmetric oper- 
ator. We show that the deficiency indices of this symmetric operator are (4,4). The 
construction is based on the observation that the coefficients if)(x) standing at singular 
terms in the asymptotic expansion in the variable Xo at a vortex of the Green function 
S z (x,x ) belong to the deficiency subspace with spectral parameter z. Here we make 
use of the explicit knowledge of the spinless two- vortex Green function 5 z (x, x ). 5 

The next and final goal which naturally follows is a construction of the two- vortex 
Green function for the Pauli Hamiltonian with the aid of the Krein's formula. Even 
this problem is solved explicitly. Surprisingly many features can be again derived from 
the asymptotic analysis near a vortex. 

The paper is organized as follows. In Section HU we summarize some basic facts 
and formulae concerning the spinless AB Hamiltonian with one and two vortices. In 
Section Mil we introduce the Pauli operator with one and two AB fluxes and derive 
the boundary conditions at a vortex defining the spin up and down components of 
the Pauli Hamiltonian. In Section HY1 we provide a basic asymptotic analysis near a 
vortex of functions from the deficiency subspaces as well as that of the spinless Green 
function. In Section El we construct a basis in the deficiency subspaces. Section ED is 
devoted to the application of the Krein's formula to our problem. 

II. Preliminaries. The AB Hamiltonian for a spinless 
particle 

The AB Hamiltonian with one vortex and describing a spinless particle, Hq, was 
introduced in Ref. [H and studied in a long series of papers by many authors. For 



2 



example, one can consult Ref. El for some mathematical details. It acts in L 2 (IR 2 , d 2 x) 
and is nothing but the selfadjoint operator associated to the closure of the positive 
quadratic form 



d Xl - i 



. 0.X2 



\x\ 



+ 



d xo + i 



\x\ 



<P 



d x . 



(1) 



defined on the space of test functions D(IR 2 \{0}). In other words, H is the Friedrichs 
extension of the corresponding symmetric operator with the domain 2)(IR 2 \ {0}). 
Owing to the gauge equivalence we can assume that a 6 (0, 1). 

We shall use the polar coordinates (r, 9) with the angle 9 G (—11,71). This implies 
a cut along the negative X\ half-axis. Sometimes it is convenient to apply the unitary 
operator 



{U a <p)(r,9) 



<p{r,9) 



and work with the unitarily equivalent operator 

H = U a H U-\ 

In particular this unitary transformation is useful when constructing the Green func- 
tion. This means that 

Dom(#) = U a (Dom(H )) . 
Formally, differential operator, 

H = -A . 

The domain of H is determined by the boundary conditions at the cut, namely 



^2tt i a 



if>(r,-ir), d r tp(r,7r) 



JI-k i a 



d r ip(r, -7r) 



(2) 



In addition, one should take care about boundary conditions at the vortex. As an- 
alyzed in Refs. OHH the domain of H is characterized by the boundary condition 
ip(0) = 0. Since ip(r,9) = exp(ia9){p(r,9) the same is true for Dom(if), namely the 
boundary condition at the vortex reads -0(0) = 0. 
The generalized eigenfunctions of H, 



J\ n+a \(kr)e 



i(n+a)Q 



fc>0,n6Z 



form a complete normalized set, 



00 ^ 

Jyikx) JJky) kdk = — S(x — y) . 

x 

This makes it possible to write down the Green function and the propagator as inte- 
grals, 



S*M;r„,9„) 



2tt ^ 



(n+a)(8-e ) 



J\n+a\ 

{kr) J\ n+a \{kr ) 
k 2 — z 



kdk 



(3) 
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and 



1 Z" 00 
% t (r,9;r ,e ) = —y2 ei{n+a){e ' e ° ] / ^ lkH J\n+a\(kr) J ]n+a] (kr ) kdk . (4) 
Z7V nez Jo 



They are related by the Laplace transform, 



Sz(r,6>;r , 9 



e xt X„ it (r,9;r ,9 )dt. 



Starting from (JH) one can derive the following formula for the propagator^ 

1 ( 1 



X t (r,9;r ,9 



1 

,27r i a 



g— 27r i a 

sin(7r a) 

7T 



4vr i t 6XP V Ait 



X — Xq\ 



(5) 



47T 



-as+ia(6»-6» ) 

_|_ g-s+i (e-o ) 



ds . 



where 



\x — x \ 2 = r 2 + r 2 , — 2r r 008(6* — 6> ), i?(s) 2 = r 2 + r 2 + 2r r cosh(s), 
and the phase factor in front of the first term depends on whether 

9 -9 e (-7T, 7r), (tt, 27r) or (-27T, -7r). 
The Laplace transformation results in a formula for the Green function, 



%(r,9;r ,9 ) 



1 

g27r i a 

g — 2-7T i Q 

sin(7ra) 

7T 



(6) 



2tt 



-zR(s] 



e -as+ia(9—9o) 
I _|_ e -s+i(6»-6» ) 



ds. 



The second term on the RHS of © can be given still another form with the aid of the 
identity 



K w {a)K-- lT {b) 



dr 



sin(7r(a + it)) 
K ^a 2 + b 2 + 2ab cosh(u) ) -j- 



-a(u—i c 



+ e - u+u 



du 



for a > 0, b > 0, < a < 1 and \<f>\ < n. This way we get 

1 



9Jr,9;r ,9 Q 



2tt 
sin(7ro;) 
2tt 2 



-z \x 



K ir (\f^zr) K_ iT (y/^zr ) 



sin(7r(o; + i r)) 



(7a) 
dr 



for - 6»o G ( — 7T, 7r), 
S 2 (r,0;r o , O ) = < 



2-7T i a 



— i^o (v 73 ^ |ac - aro|) 



(7b) 



sin 7r a) 



for 9 - 9 e (n,2ir), and 
Sz(r,0;r o ,0o) = 



2tt 2 



— 2w i a [ 



-6>o-2tt)t 



sin(7r(a + i r)) 



dr 



sm(7T o; J 
2tt 2 



|x — Sol) 



(7c) 



r ) 



-0 o +27r)T 



sin(7r(a + i r)) 



dr 



for 0-0 o e (-27r,-7r). 

Despite of this threefold description depending on the value of — 0o the Green 
function should be continuous, even real analytic, in its domain of definition if x ^ Xq. 
Checking the limits from the right and left for — O = ±7r one finds that the continuity 
is guaranteed by the identity 

/oo 
K iT (a) K_ iT (b) dr = n K (a + b) for a > 0, b > 0. 
-oo 

Let us add a remark on deficiency subspaces. First we recall a general and easy to 
verify fact. Let A be a selfadjoint extension of a symmetric operator X. Denote by 
N(z) = Ker(X* — z) the deficiency subspaces, Im z ^ 0. Then it holds 

/ e N(uO / + (z - w)(A - z)- 1 / e N(z) . 

This can be illustrated on our problem. We choose H (the one-vortex AB Hamil- 
tonian defined by the boundary conditions (J2J|) for A, and X is a restriction of H 
obtained by requiring the supports of functions from the domain of X to be separated 
from the singular point (the origin). The deficiency indices are known to be (2,2). 
For a basis in N(z) we can choose the vectors 

^_ M (r,0) =K^ a ( v ^r)e i(a -^ e , Vo,*(r, 0) = K a (V=Zr) e iaf> . (8) 

Here zGC \ R +j Re y/—z > 0. As shown in Ref. El it holds true that 



S z (r, 9; r , O ) ip ,w(r Q , O ) d0 o ) r dr 
1 



Z — W 



hence 



Vto,™ + {z-w)(H - z)~Vo 



(v 7 -^ r) - K Q (v 711 ^ r) 



(9) 
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Similarly, 



l-a 

-z 



+ (* - «0(# - *)~V-i,« = ^^z^J ^-m • ( 10 ) 

Let us now focus on the case of two vortices but still considering a spineless particle. 
The vortices are supposed to be located in the points a = (0, 0) and b — (p, 0), p > 0. 
Let (r a , 9 a) be the polar coordinates centered at the point a and (r&, #&) be the polar 
coordinates centered at the point b. The two- vortex AB Hamiltonian H is the unique 
self-adjoint operator associated to the quadratic form 

/ (\(-id Xl -AM 2 + \(-id X2 -AM) 2 d 2 x, (11a) 

where 

A = -ad9 a - (3d9 b . (lib) 

Again, owing to the gauge equivalence, we can assume that a, /3 G (0, 1). 

Also in this case one can pass to a unitarily equivalent formulation. The plane is 
cut along two half-lines, 

L a = (— oo, 0) x {0} and L b = (p, +oo) x {0}. 

The values 9 a = ±7r correspond to the two sides of the cut L a and similarly for 9 b 
and Lb. The geometrical arrangement is sketched in Fig. [TJ The unitarily equivalent 
Hamiltonian H is formally equal to —A and is determined by the boundary conditions 
along the cut, 

^(r a , 9 a = ir) = e 27ri >(r a , 9 a = -tt), <9 r ^a, 9 a = vr) = e 2wia 9 r ^(r, 9 a = -tt) , 
ij{r bl 9 b = vr) = e 2 ^(r b ,9 b = -vr), d r ^{r h ,9 h = tt) = e 27ri X^(r,^ = -vr) . (12) 

In addition, one should impose a boundary condition at the vortex, namely ip(a) = 
ip{b) = 0. 

A formula for the Green function of the Hamiltonian H is known also in the case 
of two vorticesP For a couple of points x, x G IR 2 \ (L a U L b ) we set 

Ca = lorCa = e 2 ^ Q orCa = e- 2 ^ 

depending on whether the segment XqX does not intersect L a , or ~XqT intersects L a and 
Xo lies in the lower half-plane, or XqX intersects L a and Xq lies in the upper half-plane. 
Analogously, 

0, = lor0, = e 2wi/3 or0, = e~ 2 ^ 

depending on whether the segment x^x does not intersect L b , or x^x intersects L b and 
xo lies in the upper half-plane, or xqx intersects L b and xq lies in the lower half-plane. 
Furthermore, let us set 

Ca = e ia ^, ( b = e [ ^ b where T) ai T/ b G {0, 2 7T, 2 7r}. 
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Remark. Notice that if ( a ^ 1 then necessarily £b = 1 and vice versa. 
The formula for the Green function reads 



S 2 (X,X ) = (a(b— KJ^Z \x - X \) 
ZTT 



Ca 



sm{Tra) 

2tt 2 
sin(7r j3) 

2tt 2 



- \ 1 

r 0a ) - 



^(Sa-doa-Va)r 



rot) ~ 



sin(7r(a + ir)) 



7, n>2 



sin(7r(/5 + ir)) 
1)"/ K iTn (V^r)K iiTn _^ Tn) (y^Ip) 

sin(7T(j n ) exp(0T n 



■dr 
dr 



x ... x K i{Tl „ T2) [4^~z p) K_ iTl ( 
sin(7rcr n „i) 



x 



X 



X 



'•■)) . , , . V, 

7r sm(7r(cr n + rr n )) 
sin(7rcr 2 ) sin(7ro"i) exp( 



-e Q n) 



(13) 



dV. 



7rsin(7r(cr n _i + ir n _i)) 7rsin(7r(cr 2 + ir 2 )) tt sin(7r(cr 1 + ir x )) 

Here the sum n>2 runs over all finite alternating sequences of length at least two, 
7 = (c n , c n _i, . . . , ci), such that for all j, c, G {a, 6} and Cj ^ Cj + i, and <7j = a (resp. 
/3) depending on whether Cj = a (resp. 6). In addition, (r, #) are the polar coordinates 
of the point x with respect to the center c n , (r ,9o) are the polar coordinates of the 
point Xq with respect to the center Ci(the dependence on 7 is not indicated explicitly). 



III. The Pauli Hamiltonian with AB fluxes 

According to the Aharonov-Casher ansatiP the two diagonal components of the Pauli 
Hamiltonian with the third component of spin equal to ±1/2 can be factorized, 

H ± = (p-A) 2 T B = PIP ± , 

where 

P± = (pi-A 1 )±i(p 2 -A 2 ). 

Using the complex coordinate z = X\ + i x 2 one can rewrite the Pauli Hamiltonian as 
follows: 

H + = 4 (_i^_^)(_i^_^), 
H = A{-idz-A-){-id z -A z ). 

We start our discussion from considering the situation with one vortex. Then 

a, , , . ,„ ia 91 fl l9i fl iaz ,z ia / dz dz\ 

A = — x 2 dxx - Xl dx 2 ) = -add = — e~ 2ie de 2l(> = — - d- = — — . 

r z 2 2 z z 2 \ z z J 

Hence 

ia ia 

Az = Yz' A ~ z = ~^ 
7 



and we can write 

In fact, these are formal expressions. More precisely, the operators are defined as the 
unique selfadjoint operators associated respectively to the positive quadratic forms 



q+O) = 4 



\ 2 

d 2 x, q_(^) = 4 



d 2 x. 



(14) 



with their natural maximal domains of definition. 

Since the magnetic field vanishes on M 2 \ {0} the operators coincide with the 
spinless AB Hamiltonian H on the domain 2)(IR 2 \ {0}) (2) is the space of test func- 
tions). This means that all three operators H + , H~ and H are selfadjoint extensions 
of the same symmetric operator X. From Refs. ElandEJit is known that all selfadjoint 
extensions can be described by appropriate boundary conditions at the origin. The 
method used to derive the boundary conditions was inspired by the description of 
point interactions in the plane given in Ref. Let us also note that analogous bound- 
ary conditions have been derived in Ref. [HI for the model with additional homogeneous 
magnetic field while the case of Dirac-Weyl operator is discussed in Ref. El 

To describe the boundary conditions one introduces four functionals, 



2k 



l imr i-«_ / (p (r,e)e i9 de, 
rio 2ir J 



lim r 



-l+a 



27T 
2tt 



2 Ti 



<p(r,6)e l6 d6 - r -1 ^^ 1 ^) 



o 



limr°— / <p(r,8)d6, 
no 2n J 

f2w 



\imr- a (—[ <p(r,6)d6 -r- ® ^)) . 

rlO \27l J J 



Each boundary condition is determined by a couple of matrices Ai,A 2 G Mat(2, C) 
fulfilling (the symbol (Ai,A 2 ) designates a 2 x 4 matrix) 



rank(Ai, A 2 ) = 2, A 1 D~ L A 2 : = A 2 D- l A\ 



where 



D 



1 - a 
a 



The boundary condition takes the form 



A, 



A, 



Two couples of matrices, {Ai, A 2 } and {A[, A 2 }, determine the same boundary con- 
dition if and only if there exists a regular matrix G G GL(2, C) such that (A[, A' 2 ) = 
G(A 1 ,A 2 ). 
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For example, the domain of the spinless AB Hamiltonian Hq is determined by 
the boundary conditions at the vortex ^ 1 (lp) = = and so by the couple of 

matrices 

(I o\ . fo o N 



\) ' z \0 

We wish to derive the boundary conditions for the Hamiltonians H + and H~ . Ac- 
cording to the well-known construction, the operator A associated to a semi-bounded 
quadratic form q is determined by the condition 

V/ G Dom(A) C Dom(q), G Dom(q), (<p, Af) = q(^, /) . 

This is to say that / G Dom(q) belongs to Dom(yl) if and only if there exists g G "K 
such that the equality (</?, g) = q(ip, f) holds true for all ip G Dom(q). In that case g 
is unique and Af = g. We are going to apply this prescription to the quadratic forms 
(|T3|). This amounts to integration by parts. 

More precisely, the Green formula implies that 

/ (d z f)gd 2 x = -[ f(d z g)d 2 x-\im^ [ 2 \fg)(acos(9),asm(9))e- ie d9. 
Jr 2 Jr 2 a ±° z Jo 

Thus one finds that / G Dom(X*) belongs to Dom(H + ) if and only if for all <p G 
Dom(q + ), 

lima / (jp( dz- ^z)f)(a cos(#),a sin(#)) e" ie d# = 0, 

or, when expressing (z,z) in the polar coordinates, 

lim / (lp(rd r +id e -a)f)(acos(6),asm(6))d6 = 0. 

Any / G Dom(X*) asymptotically behaves like 
/ = ($^(/) r~ 1+a + $2 *(/) r l ~ a ) e~ ie + ($?(/) r" a + $°(J) r a ) + regular part. 
Hence 

(r <9 r + i d e - a)f ~ 2(1 - a)$ 2 *(/) r^e" 16 - 2a $?(/) r~ Q + . . . . 
Notice that 

(r 9 r + i d e - a)r~ 1+a e-' ld = {rd r + id e - a)r a = 

and so any function of the form r~ 1+a r](r)e~ l6 or r a 7](r), with 77 G C°°(IR + ), ry(r) = 1 
in a neighborhood of and r/(r) = in a neighborhood of +00, belongs to Dom(q + ). 
Therefore a sufficient and necessary condition for / to belong to Dom(H + ) is 

^\f) = *?(/) = 0- (15) 
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The corresponding couple of matrices can be chosen as 



• o i\ . /o o 



oy ' z ^i o, 

The other component of the Pauli Hamiltonian, H~ , can be treated similarly. One 
finds that / G Dom(X*) belongs to Dom(H~) if and only if for all <p G Dom(q_), 

lim / (jp (r d r — id s + a)f)(acos(8),asin(8)) d9 = 
which turns out to be equivalent to 

$r 1 (/) = ^(/) = o. (i6) 

The corresponding couple of matrices can be chosen as 



0J ' z \0 1 

The generalization to the case of several vortices is quite straightforward. One 
simply imposes the above derived boundary conditions at each vortex. Let us consider 
the case of two vortices. For the sake of simplicity we assume that the vortices are 
a = (0, 0) and b = (1, 0). The Pauli Hamiltonian formally reads 



2\z z-lJ ) \ 2\z 2-1 



2\z z — 1/ J \ 2\z z — 

We still assume that < a,/3 < 1 (in virtue of the gauge equivalence). 

The Pauli Hamiltonian with two vortices is known to have zero modest They can 
be computed with the aid of the Aharonov-Casher ansatz since it effectively enables to 
reduce the second order differential equation to a first order one. Explicit solutions are 
even known in some essentially more complicated situations (see for example Ref . ITTJ) . 
Just for the sake of illustration let us verify that the zero modes actually satisfy the 
above derived boundary conditions (JT5|) or (jUJ). 

If a + (3 < 1 then the function 



<p(z) 



1 



i/3 



z(z 



is L 2 integrable and solves 



*-i(2 + i | T ,|„-o. 



So it is a zero mode of H + . It is elementary to compute its asymptotic behavior for 

<P = r- 1+a e- ie ° + (l - - §Ce- 2ie ° + O(r^). 
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Hence <p obeys (jT5|) . The boundary condition at the vortex b is analogous. 
Similarly, if a + (3 > 1 then 

\z\ \z — 1| 

is a zero mode of if - and 
Hence <p obeys (JTHJ). 



IV. Asymptotic behavior near a vortex 

Our first task in this section is the asymptotic analysis of functions from a deficiency 
subspace. To simplify the discussion we shall use the symbol 0(r 7 ) in a sense somewhat 
weaker than it is common. The equality f(r,9) = 0(r 7 ) for r j will mean that 
f(r, 9) = J2nez fn(r)e ine and for all n it holds f n (r) = 0(r 7 ). 

Lemma 1. Assume that R > 0, z G C \ R + , < a < 1 and cp G L 2 (B(0,R), d 2 x) 
satisfies in the weak sense the differential equation 

{Y-z)<p = 

on B(0,R) \ {0} (the disk centered at with the radius equal to R) where (using the 
polar coordinates (r, 9) ) 

Y = - e - ia6 Ae iad = -(d Xl + iad Xl 9) 2 - (d X2 + iad X2 9) 2 
d 2 Id 1/0 \ 2 
dr 2 r dr r 2 \ 39 

Then there exist constants cq, do, c_i, d-\, such that 

p(r,9) = c r- a + d r a + (c- ir - 1+a + d_ 1 r 1 - a )e- ie + 0^) for r 10, (17) 
where 7 = min{2 — a, 1 + a}. 

Proof. For all n G Z, 77 G 2)((0, i?)) (the space of test functions) it holds true that 

= ((Y - z),(r)e'«<» = -( (|j + i| + z - ^^-)n(r) 
Hence 

00 

<p{r,0)= fn(r)e ine , 
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where 

in the weak sense. This implies that the generalized derivative d r (r d r f n (r)) belongs 
to L\ QC ((0, R)) and consequently /„ G AC 2 ((e, R)) for all < s < R. Therefore 
necessarily f n ( r ) is a linear combination of the modified Bessel functions, 

fn(r) = a n K n+a (y/^z r) + b n I\„ +a \(V-z r) . 

Let us recall the asymptotic behavior of the Bessel functions. If < v and v ^ N then 

and 

A -„ (r) .£M (r(1 + 0(r2)) _£(Lzii ( i r(1 + 0M) . 

This implies that f n G L 2 ((0, R), r dr) if and only if either a n = or n G {0,-1}. 
This is to say that a n can be nonzero only for n — 0, —1. So if n 7^ 0, —1 then 
/ n (r) = 0(r' n+a '). This proves the lemma. □ 

Let id be the two- vortex spinless AB Hamiltonian defined by boundary conditions 
(|T2l . The symbol X below stands for the symmetric operator obtained by restricting 
the domain of H so that functions from Dom X vanish in some neighborhood of the 
vortices. The deficiency subspaces are denoted by = Ker(X* — z). 

Corollary 2. // z G C \ R+ and if) G 3\T(z) £nen taere exist constants c at o, c a -\, c^o, 
Cb-i, such that 

ij(x)=c afi r; a e la ^+c a ^ 1 r~ 1+a e l(a ~ 1 ^+o(l) for r a j 0, (18) 

and 

0(x) = c 6 , e i/3 ^ + c b ,_! r" 1+/3 e^" 1 ^ 6 + o(l) for r b | 0. (19) 

Proof. The property V € N(z) means that V £ L 2 (B(0,R), d 2 x), (-A - 2;)^ = 
on M 2 \ (L a U L b ) in the weak sense and ip satisfies the boundary conditions (fT2|) on 
L a U L b . Then the function exp(— ia9 a )ip obeys the assumptions of Lemma [T] and 
relation (fTTj) implies (fT8|) . Relation (fT9|) can be shown similarly. □ 

Corollary 3. Assume that z G C \ R +3 ^ G N(z) and = V'(^) = 0- Taen 

■0 G Dom(dd) and hence ip = 0. 

Proof. We use once more the fact that exp(— i a 9 a )ijj obeys the assumptions of Lemma[T] 
and hence 

exp(-ia6 a )^(x) = c r- a + d r« + {c^ r - 1+a + d^r^e^ + 0(r^) for r a [ . 

Since i])(a) = it holds c_i = c = 0. Let {7 be the unitary operator on L 2 (M. 2 , d 2 x) 
acting via multiplication by the phase factor exp(ia9 a + i(36 b ). Then ip = exp(— ia6 a — 
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ij39b)ip belongs to Ker(X* — z) where X = U 1 XU. The function 6b(x) is real analytic 
in a neighborhood of a and 

9 b (x) = sm(9 a ) - + 0(r 2 a ) for r a j . 
P 

A straightforward computation gives the asymptotic behavior of (p and one finds that 



^\<p) = C-l, $2%) = d-1 + $?(</?) = C , $»M = d, 



1 



2p' ~ lvr/ ~ u ' ~ zvr/ ~° 2p ' 

So one finds that the boundary condition §i l (ip) = $i(<p) = is satisfied at the 
vortex a. Analogously, the same boundary condition is fulfilled at the vortex b. As 
recalled in Section ITTT| these boundary conditions determine the domain of H . Hence 
(p G Dom(ifo) and ip G Dom(iJ). But H is positive, z G" M+, and therefore Dom(if ) fl 
Ji(z) = {0}. This shows that ip = 0. □ 

Further we are interested in the asymptotic behavior near a vortex of the Green 
functions and (fT3|) . It is easy to see that in the spinless case the Green function 
vanishes in each vortex. For example in the case of two vortices it holds true that 
Sz(a,a;o) = 0. This can be derived from ffTTTf) with the aid of the relation 

K iT (r) -> 7i 5 (t) forr|0 (20) 

and some simple combinatorics. It is also obvious that 

K (y/^z \x - x \) = K (y/^zr 0a ) + 0{r a ) for r a | 0, 

(here r 0a = \a — x \) and 

K iv (y/^r b ) =K iu (y/^p) +0(r a ) for r a j 0. 

To get an additional information we shall need an asymptotic formula for the 
integral 

f°° ir I l — v (l — ^ sin(7ra) exp(9 a r) 

/ K iT (V-zr a ) K i{v _ T) (y/-zp) —— — —--dr. 21 

J-oo 7r sm(7r(a + ir)) 

Such an asymptotic analysis can be carried on with the aid of the following lemma. 
Lemma 4. Suppose that r > 0, \9\ < n and < a < 1. Then 

i + e -u+s+w sinfvra) 1 - a \2/ 



T(l - a) (t 



a \2 



Q a e~«^ + Z(r,u), (22a) 



where 



Vr G (0, 1), Vm G M, \Z(r,u)\ < Krcosh(u) (22b) 
and K depends on 9 and a but does not depend on r and u. 



13 



Proof. The LHS of (l22al> equals 



n a)(s— poo p -a(s+u— i 



p -rcosh(s) ^^ ds+ / e -rcosh( S ) ds _ 



Therefore it suffices to study integrals of the form 

e -rcosh(s)_f ^ _ e y(u-i6) / p -r cosh{s)-~/s ( j g 

1 + e - s +"- ie y u 

00 -(7+l)(s-M+i6») 

e -r cosh(s) ds / 24 ^ 

]_ _|_ g— s+it— 10 v ' 

for < 7 < 1. The second integral on the RHS of (j2*IJl can be treated easily and one 
finds that 



e -r cosh (s) e ds = e / e ds + Zi(r, u) 



where Zi(r, u) satisfies estimate fj22b(l . To treat the first integral on the RHS of J2 
we note that 

CO ^ , 
-SIX i \ 



fc=0 

and therefore 



/"OO 1 

/ e" r cosh («)-7« ds = V - 2^~ V +fc r (-7 - k, -) (- ~ 
J k=Q k. V 2/ V 2 



T(l - 7) /r\7 1 7 „ 
'- +- + - 9 r + 0(r 2 ). 



7 \2/ 7 1 — 7 2 

Furthermore, 

/ e -^osh(,)- 7 , ds = + z , s 

Ju 7 
where Z 2 (r,u) satisfies estimate (|22b|) . Thus we have derived that 

J u i + c -s+ u - 1 e 7 V2/ 7 l + e- s - ie v ' 



(25) 

where Z 3 (r, w) satisfies estimate (|22bj) . To conclude the proof it suffices to apply 
to the both integrals in and to take into account that 

00 e _(l_ a )(s+i6») roc p - a ( s -i6») roc - a ( s -W) n 



ds + - — ds = — ds 



/ 1 + e- s - ie ' J 1 + e- s+w J_ 00 1 + e~ s+ie sin(Tra) 

for 101 < 7T. □ 
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Corollary 5. Under the same assumptions as in Lemma^ it holds true that 



K iT {V-zr) A i(v _ T )(V -z p) — dr 



7T 



sin(7r(a; + it)) 



K w {yT^zp) 



sin(7ra) T(a) 



7r 1 — a 
siniir a) r(l — a) ( sf—z r 



-z r 



l-a 



A^ K . wl+a ^ zp ) 



71 



i a9 K iu+a (^p)+0(r) 



(26) 



for r I 0. 
Proof. Using 



K iT (a) = - I e iST - acosh(s) ds for a > 0, r G 



(27) 



and applying the equality 



ciT fu-,) exp(gr) ^ = 2 e -°(»— *> 
sin(7r(a + ir)) 1 + e - u+s+w 



we find that (|2*T1) equals 



sm(7raj 
2tt 



-2 p cosh(tt)— i v u 



-z r cosh(s) 



-a(u—s—i0) 



1 + e 



-U+S+10 



ds du 



Now it suffices to apply (|22j) to the inner bracket and then to use the integral form 
(1271) in the reversed sense. □ 



First let us apply (|26|) to the case of one vortex. In fact, the following observa- 
tion about the asymptotic expansion of the Green function near the vortex will 
be crucial for the subsequent analysis. We get either the asymptotic expansion of 
S z (r, 9; r , 6 ) for r j 0, or, since in general it holds true that 



9*(r, 8; r , O ) = Sz(r , ] r, 6), 
the expansion for r j as well, namely 



Sz(r,9;r ,6 ) 



sin(7ro;) r(a) 
2tt 2 l-a 



l-a 



K-i+c^x/^r) e 



i(a-l)(0-0 o ) 



(28) 



(29) 



+ 



sm(n a) T (1 — a) fy/—zro 



2tt 2 



a 



K a (V=z~r) e ia{d - eo) + O(r ). 



One observes that the coefficients standing at r$ a e l ( a l ' e ° and r$e ia9 ° are respec- 
tively proportional to 

K„ 1+a {y/=z~r) e i{a - 1)e and K a (^r) e ia6 '. 
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But these functions are nothing but the basis functions in the corresponding deficiency 
subspace, see (JEJ). 

Next we shall consider the case of two vortices. Applying (|2~K|) to (TTTT|) we get 



sin(7ro;) T(a) 



2tt 2 I- a 



-zr r , 



l-a 



+ 



sin(7r a) T(l — a) 



2tt 2 



a 



-z r, 



L a (x ) + 0(r a ) (30a) 



for r a [ where 

L v {x Q ) = K u (V^r 0a ) e-- e - + V (-1)"" 1 / K iTn _ 1+I/ (^p) 

7,n>2, Cn =a • /Kn - 1 

x ^(^-^(v^p) x ... x K i(ri _ T2) (V^p) ^^(v^iro) 

sin(7rcr n _i) sin(7ro- 2 ) sin(7rai) exp(-0 O Ti) x 



x 



7rsin(7r(a n _i + ir n _i)) 



X ... X 



7rsin(7r(cr 2 + ir 2 )) 7rsin(7r(cri + in)) 



(30b) 



(and, again, (r , 6* ) are the polar coordinates of the point x with respect to the center 
Ci). The convergence of the series in ()30b|) will be discussed later in Section Ivl 



V. Deficiency subspaces for the case of two vortices 

In this section we are going to construct a basis in the deficiency subspaces in the 
two-vortex case. So H designates the two-vortex spinless AB Hamiltonian described 
by the boundary conditions (fT2|) . X is the symmetric operator obtained by restricting 
the domain of H as described in Section HXl and 'N(z) = Ker(X* — z) is a deficiency 
subspace. 

Asymptotic expansion (|29|) for the one vortex case suggests that also in the two 
vortex case one may extract from the Green function a basis in the deficiency subspace. 
From (pffij) and (|28|) on derives immediately a candidate for such a basis. It is formed 
by the functions 

oo 

ipu,v,z( x ) = ^ S n {u, u, z; x) , (31a) 

n=0 

where 

S (u,u,z;x) =K v (V=z~r u ) j u6 \ (31b) 

S n (u, is, z; x) = (-1)" f K iTn (V=i r n ) (31c) 

x K i(Tn _^ Tn) {^z~p) x ... x Kifa-Tt^V^zp) K„ iTl - v (^zp) 
sin(7ro- ra ) exp(6 l w r n ) sin(7rcr n _i) sin(7rcri) 
7rsin(7r((T n + ir n )) 7rsin(7r(cr n _i + ir n _i)) 7rsin(7r(o-i + in)) 
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for n > 1, the indices are restricted to the range 

u G {a, b}, v G {uj — 1, uj} where to = a if u = a, and to — (3 if it = b, (31d) 

and to each n E N one relates the unique alternating sequence (c„, c n _i, . . . , ci), Cj G 
{a, 6} and Cj ^ Cj + \, such that c\ ^ u. Correspondingly, cjj = a if Cj = a and <7j = (3 
if Cj = b. As usual, (r n ,9 n ) = (r Cn ,8 Cn ) are the polar coordinates with respect to the 
center c n , (r c ,9 c ) are the polar coordinates centered at the point c. 

Let us show that the series (|31a|) actually converges. In the Hilbert space L 2 (IR, dr) 
we introduce the vectors 

fu,z( x '-> T ) = K iT (^f^zr u ) exp(6> M r) . sm ( 7rcr ) . Q v , z {r) = K^ lT A\f^zp) , 

v ' 7r sm(7r(cr + it ) ) v ' 

and the operators & z and D u with the generalized kernels 

_ , \ T s ( I \ ^ / \ sin(7ro") . 

£ z (/x, uj) = \y/-zp) , w) = —7 — — -tt o(/x - w), 

y v 7 7r sm(7r(cr + i/i)J 

where 

u G {a, 6}, a = a if w = a, and a — (3 if u = b. 

For k G {a, 6} let t> be the complementary vortex, i.e., {it, t>} = {a, b}. For the sake 
of brevity we shall use the matrix-like notation in the following paragraph. Thus the 
transposition will in fact indicate an integration , i.e., f T Q = L f(r) g(r) dr. 

We can rewrite the summands in equation (|31a|) using this notation (here n > 1), 

San_i(u,i/,z;a:) = -f^(x) T (^ S) u ^D„) n_1 0^ , 
S 2n (w,z/,2:;x) = f^^^D^^D^^D^)""^^. 

These formulae make it possible to estimate the summands. Note that & z acts as a 
convolution operator and so it is diagonalized by the Fourier transform. Since 

e iXT K iT (a) dr = ^e'^H^) 



IT 

OO 



we get 

|| jy = Tie- Re{ ^ )p . 

The operator D u is already diagonal. Therefore 

sin(7ro") 



\® u \\ = SUp 



7r sin(7r(cr + i/x)) 



1 

7T 



Jointly this implies that 

IS^fai/,*;*)! < ||M^)IIII^He- ( ^ 2)Ha(V=i)p , 
\S 2n (u,u,z;x)\ < Hf^^llll^lle-^- 1 ^ 6 ^". 
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The estimates show that the series (|31a,j) converges absolutely at least as fast as a 
geometric series. Even one can rewrite the formula for ip u ,u,z(x) in a compact form, 
namely 

^u,u,z{x) = K v (^r u )e lv9 * (32) 

vl av.z- 

Here the inverse operator (I — & Z 'D U & Z '2) V )'~ 1 exists with the norm estimated from 
above by (1 — exp(— 2 Re(-\/— ~z) p)) _1 . 

Altogether we get four functions: ip a ,a-i,z, i>a,a,z, ipb,0-i,z an d ipb,p,z- Our goal is to 
show that they actually form a basis in the deficiency subspace. Obviously 

(A + z)Vv,* = 



since 



{A + z)K v (^r)e ±iv6 = on R 2 \ {0} 

for all v E C, z E C \ R + , and therefore all the summands satisfy the equation 
(A + z)S n (u, u,z)=0 in the domain M 2 \ (L a U L b ). 

Let us verify that ip u ,u,z obeys the boundary conditions (fT2|) . For the sake of 
definiteness we shall consider the function ip a ,u,z% v E {a — l,a>}. Firstly we shall show 
that 

e- ia ^ a ,,,z\ 6a=n -e ia ^ a , u , z \ ea= _ n = 0. (33) 

If n = 2m — 1 is odd then c n = b. Moreover, if 6 a = ±n then 9 b = and r b = r a + p. 
Hence 

e~ ia7T S 2m -i(a,v,z)\ ff -e ia7r S 2m -i(a, v,z)\ a _ w 

= / ^ir 2m _! (V^z(r a + p)) x . . . x K iiTl _ n) (V^zp) K- iTl - v {y/=z p) 

x 2isin(7ra)sin(7ro- 2m -i) x x sin(Tro-i) d 2m-i r 
7rsin(7r(o- 2m -i + ir 2m -i)) 7r sin(7r(<7i + in)) 

If n = 2m, m > 1, is even then c n = a and 

e" ia7r exp(7rr n ) - e 1Q7r exp(-7rr n ) = -2i sm(n(a n + ir n )) 

hence 

e~ ia7T S 2m (a, v,z)\ _^ - e ia7T S 2m (a,u,z)\ _ 

K iT 2m {V^r a ) ^i(r 2m _ 1 -r 2m )(v /Z i P) X ... X ^^(V^i^j) 



x 2isin(7ra) sin(7ro- 2m -i) ^ ^ sin(7r n) d2m ^ 

7T 7rsin(7r(o- 2m _i + ir 2m _i)) "' 7rsin(7r(a 1 + in)) 

The integration in r 2m can be carried on with the aid of the identity 

K iT (a) K i[v _ T) (b) dr = vr K iv (a + b) for a > 0, b > 0. (34) 
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This way we get the equality 

e- iaw 5 2m (a,z/, z)\ g _ 

i u a — ■ ■ 

= -(e" iM Vi(a l v,z)\ 
valid for all m > 1. Obviously, 



e la7r S 2m (a,^z) 



3 ia *SW-i(a,ML a= _J 



e- ia "S o (a,v,z)\ 0a=n -e ia *S (a,v,z)\ 

The last two equalities imply (f33|) . 
Similarly one can show that 

e~ i/3?r w> I _ e i/3^w, I 

I (7Jj — /( 



0. 



y b =-ir 



0. 



(35) 



Equality (|3~i|) again turns out to be useful but this time when treating the odd sum- 
mands. With its aid the dimension of the integration domain is reduced by 1 and one 
obtains the equality 



e lf3n S 2m -i(a, v, z)\ - e l/3n S 2m -i{a, u, z)\ 
= -(e- i/3 ^ 2m _ 2 (a, u, z)\ - e i ^5' 2m _ 2 (a ) u, z)\ ) 



V b = -W 



valid for all m > 1. This shows (pTHIi . 

Finally we note that the remaining two boundary conditions, 



<9r n 



dr h 



<9r a 
dr b 



0. 



0. 



1b=-n 



can be verified in exactly the same way. 

Next we wish to examine the asymptotic behavior of the functions ip UjI/>z near the 
singular points a and 6. We shall again focus on the functions ip a ,v,z, the functions 
i>b uz can be treated similarly. First notice that 



(36) 



Actually, for the even summands in (|31a,j) the limit x — > b just means setting r a = p. 
To treat the odd summands one applies the limit procedure (f20|) for r b — ► and finds 
that 

S2m-i(a, v, z; b) = -S 2m ^ 2 (a, u, z; b) for all TO > 1. 

This shows (JBEJ). 

Let us make this result more precise. The even summands in (|31a,|) simply satisfy 



S 2m {a, v, z; x) = S 2m {a, u, z; b) + 0{r b ) for x -> b. 
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Asymptotic behavior of the odd summands can be obtained with the aid of relation 
(j2EJ). We get (here S 2m -i{a, u, z\ b) = -S 2m -2{a, u, z; b)) 



S2mr-i(a, u, z\ x) = S 2m -i(a, u, z; b) + 



sin(7r/3) T((3) ( \f-ZT}, 



1-/3 



3 iC/3-i)e 6 



x 



X 



7i 1-/3V 2 
i +y9 (V— x ... x K- l{T1 _ T2) (yf^z p) K_ iT1 _ u (y/^zp) 



Sm(7T0 2m -2, 



X . . . X 



Sin(7T (Ti . 



7r sin(7r(cr2 m -2 + ir 2m -2)) tt sin(7r(cri + in)) 



d 2m-2 T 



+ 



sin(7r/3)r(l f^f=z~r b 



71 



(3 



j-Pt 



x 



X 



iT2m-2+/3 



(V^p) x ... x Ki^-^) ^-^-^(V^p) 



sm( 7T cr 2m _2 y 



7rsin(7r((T2 m _2 + ir 2m _2)) 
for x — > 6. Jointly this means that 



x . . . x 



sm(7r 0\ 



7rsin(7r(a"i + iri)) 



Tpa,v,z( X ) = ^2 

for x — > b where 



sin(7r \fi\) T(l - \p\) ( y/^zr b 



7T 



I Ml 



e i ^§^(a,(3-z) + 0(r b ) (37) 



oo „ 

S W) „(a,/9; z) = K^ v (-^zp) + / ^^(v^p) 



x ^i(r 2m _ 1 -r 2m )(v / ^p) x ... x ^^.^(V^zp) ^^^(v^p) (3f 
sin(7ro- 2m ) . . sin(7rcri) l2m 



x 



x . . . x 



7rsin(7r(a 2m + ir 2m )) 7rsin(7r(cr 1 + iri)) 



with (a 2m , . . . , cr 2 , Oi ) = (a, . . . , a, /3). 

The function ip a>V)Z has a singularity at the point a. Nevertheless it holds true that 

oo 

S n (a, v, z; a) = 0. (39) 

n=l 

The verification is similar to that of equality (|36|) . This time it holds true that 

5 , 2m(o, ^, ^; a) = — 5 , 2m _i(a, z/, z; a) for all m > 1. 
This shows (|39|) . A more precise result can be derived as follows. Note that 
S 2m -i(a, u, z; x) = 5 2m _i(a, u, z; a) + 0(r a ) for x -> a. 
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Asymptotic behavior of the even summands can be obtained with the aid of relation 
IjSEjl . We get 



S 2m (a, z; x) = S 2m (a, u, z\ a) - 



sin(7ra) T(a) f y/—zr a 
n 1 — a 



l-a 



X 



/ K iT2m _^ 1+a (y/^I p) X ... X K i(Tl . T2) (y/^ p) K- in - v (y/^Z p) 



sm(7rcr 2m -i 



X . . . X 



sm(7r<7ij 



7T 



sin(7r(a 2m _i + ir 2m _i)) " ' vr sin(7r(o"i + in)) 



d 2m-l 



sin(7ra) T(a) 
7T 1 — a 



—ZTn 



l-a 



X 



sin(7T(Ti) 



sml7rcr 2m -i 



X . . . X 



7Tsin(7r(a 2m _i + ir 2m _i)) 7rsin(7r((Ji + in)) 



d 2m -V + 0(r a ) 



for x — ► a. The asymptotic behavior of the Macdonald function is given by the 
formulcpl 



7T 



KJx) 



2 V2 

Finally we arrive at the expansion 

-M 



x 



2sin(v7r) VX(1 -v) Y{\ + v 



+ 0{x~ v+2 ) 



' ' + Ofc ) for < v < 1. 



2u 



r(H 



- £ 



Sill(7r |/i|) r(l - \p\) ( ^J^ZTa 



71 



\p\ 



(40) 



^ z) + 0(r a ) 



for x — > a where 

^lAjyVifi^l Pi Z) ■ 



71 



2 sinfvra 



oo „ 

/T|j2m — 1 
m=1 </K 



iT2m-l+^ 



(v^p) 



X ^i(r 2m _ 2 -r 2m _ 1 ) (v^p) X ... X K i(ri _ T2 ) (v^p) K-^-^V^Z p) (41) 

sin(7rcx 2m _i) sin(7rcri) 



X ... X 



7rsin(7r(cr 2m _i + ir 2m _i)) 

with (a 2m _i, . . . , er2,cri ) = (/?, • • • 

Remark. As a consequence one can show that 



7rsin(7r((Xi + iri)) 



d 2m -V 



^2 S n(u, V, Z; X) = ^ S n{u, ~V, Z\ x) . 



(42) 



n=l 



n=l 
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Actually, a short inspection of the above derivation shows that the functions 

oo 

71=0 

also satisfy the boundary conditions (fT2|) and solve the equation (A + z)ip u ,v,z = 0. 
Therefore the function 

oo oo 

f(x) = 22 S n(u, v, z\x)-y j S n (u, -u, z\ x) 

n=l n=l 

satisfies the boundary conditions (fT2|) as well and solves (A + z)f = 0. In addition, 
f(a) = f{b) = 0. Consequently, / G Dom(if ) and (H — z)f = 0. Necessarily / = 0. 

Lemma 6. dim N(z) < 4. 

Proof. In virtue of Corollary for any five-tuple of functions from there exists 

a nontrivial linear combination of these functions vanishing both at a and b. By 
Corollary O such a linear combination equals 0. 

□ 

Proposition 7. dim!N(;z) = 4. 

Proof. Owing to Lemma El it suffices to show that dimN(0) > 4. But in relation 
(j3T| we have constructed four functions ip a ,a-i,z, tpa,a,z, i/>b,p-i,z an d ^b,p,z from the 
deficiency subspace J{(z). The asymptotic expansions (pTfjl and (jinjl show that these 
functions are actually linearly independent. □ 

We conclude that the functions {tp a ,a-i,z, ipa,a,z, ipb,p-i,z, i*b,p,z} form a basis in 'N(z). 

Remark. Formula (|32j) is well suited for numerical computations. To give the reader 
an idea about the behavior of ip u , U)Z we have plotted l^a-i.il in Fig. [2] and i n 
Fig. El with a = 1/3, f3 = 2/3 and p — 1. Note that the former function vanishes in 
the vortex b while the latter one vanishes in the vortex a. 



VI. The Krein's formula 

We would like to emphasize once more that we are using two unitarily equivalent 
formulations. The operators H ± , H Q are respectively associated to the quadratic 
forms (Tl11) and (jTTl) . Let U be the unitary operator in L 2 (R 2 , d 2 x) acting as Uip = 
exp(ia9 a + i/39b)ip . The Green function (|T3|) corresponds to the operator H = UHqU^ 1 
defined by the boundary conditions on the cut ([T2|. 
Set 

fu,u,z = V-z 1p u ,v,z- (43) 
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Let us enumeratethe basis {f a , a -i,z, fa,a,z, fb,p-i,z, fb,p,z} in N(z) as f z , /f , / 2 4 } (in 
this order). Set p z = U' 1 ^, R z = (H - z)' 1 , Rf = (H* - z)' x . According to the 
Krein's formula 

Rf-R z = J2(M±y> k P z (fl-) (44) 

or, in terms of Green functions, 

Sf(x,x ) = S z (x,x ) + Y,(M±V' k tt(x)m^) } (45) 

where Mf is a holomorphic matrix- valued function defined on C \ M. 

An operator-valued function Rf constructed this way will be the resolvent of a 
selfadjoint operator if and only if it satisfies (Ref. UHl Chp. 5.2) 

V2GC\1, (Rf)* = Rf (46) 

and (the Hilbert identity) 

Vz, weC\8, Rf-Rf = (z- w)RfRf (47) 

(it follows from ijHJ) that Ker Rf = {0} for all z G C \ R). Let us analyze conditions 
(|46|) and (|47|) . It is straightforward to see that (f46|) is satisfied if and only if 



M* = Mg. (48) 

In equality (f52|) below we shall show that 

Vz, w e C \ R, Vj, fi + (z- w)RJi = Jl 

With the aid of this identity it is just an easy computation to show that (JITJ) is 
equivalent to the condition 

Vz,weC\R, M z - M w = (z - w) M z P(z, w) M w , (49) 

where P(z,w) is the 4x4 matrix of scalar products, 

p{z,wY k = uif k w ). 

Equality (JIHJ) was presented in Ref. HU and was applied to problems similar to ours 
for example in Refs. Upland l3l 

According to formula (|3Hj) and definition (pTT| of ip u ,u,z(%) we have 

sin(7r«) T(a) ( \/—zr 0a s ' 



a 



sin(7ra)r(l-a) ( V^r 0a \ _ iQdQa 



+ 2 ; 2 a I ) e 0a + °( r o«) (50) 
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for r 0a | 0. Using this asymptotic behavior and the Hilbert identity written in terms 
of Green functions, 

(z-w) / 9 z (x,y) Sw(y,xo) d 2 y = S z (x,x ) - 9w(x,x ) , 

we obtain an equality valid for u = a, namely 



[x,y)ip u>u , w (y) d 2 y 



(51) 



M 



This means that 



^,i/ )lu + (2;-«;)(-H"-2:) Vu,i 



for v E {a — 1, a} and u = a. The same argument naturally applies also to the vortex 
u = b. Using notation ijlHJ) we find that 



fu,v,w W^)(H z) fu,v,w fu,u,; 



(52) 



holds true for all u>, z E C \ R. 

We wish to compute the 4x4 matrix of scalar products P(z,w). Using (|28|) and 
applying the asymptotic behavior ijHUI) once more, this time to equality (|5TJ) . we find 
that the integral 



(y) d 2/ 



equals the coefficient standing at 



sin(7r |/i|) r(l - |/i|) /r B \H , ga 



2tt 2 V 2 

when taking the asymptotic expansion of the expression 



and 



*-^(v^) H (v^) H 
for a; — * t> , i.e., r v { 0. In virtue of pTil) and fpT7|) we get 

(y) d y 

1 



yj2 



-2vr 



Z — W \ \ y/—W 



= \H / / 

T M>i ,(a,/?;z) - 



lpa,»,z(y) ^b,u,w(y) d y 

= 2. 1 



z — w \ y^—w 



m 



- : \ „ I \J-w 



S Vtlt (a, P; z) 



S w/i (a,/3; w) 
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In particular, 



2tt 



\va{z) 



Im I 



—z 



M 



7 u>v (a, P; z) 



This means that, when passing to functions {/|} instead of {^^z}, 



— w) P(;z, w) 



-2tt 



-lT Q _i, a (a,ftz) 
"~ S a _i ij9 _i(/3,a;z) 
§ Q _i i/3 (/3,a; 2) 



-«T a a -i(a!,/3; z) 



(V^z) 2a 7 a , a {a,[3; z) 



V-(^f) 1 -"^' 



-\l+a-p. 



l a ,p-i((3,a;z) 
(y/^£) a+ S a>f3 (P,a;z) 



— \ l+a-0 



(v/^i) 7p-i,p-x{fi, a;z 



Sp-i,a(a,(3;z) 

,/9-l(A« 



- ( v 7 ^) " +/3 S/3,Q (a, /3; 2) 

—\Wc 



(V^z) 3> >/9 (/3, at; «) J 



(53) 



The Green functions Sf(x,x ) should satisfy the corresponding boundary condi- 
tions in each variable x, x . Let us first consider the case of H + . Recall that the 
boundary conditions which determine the domain of H + are = $o = (see (JT5|). 
Let us check the asymptotic behavior of Sf(x, xq) for x — * a. Asymptotic behavior of 
3z{x, Xq) is given in ([571 and asymptotic behavior of /j(^o) follows from pH|l and (pTfj) 
jointly with definition l|4*3*j) . The condition = means that the coefficient standing 
at ( r oa/2)~ a exp(— ia^oa) vanishes. This term occurs only in the asymptotic expansion 
of fl(xo) and so 

j 

The set of functions {/|} is linearly independent and thus we get a condition on the 
matrix M+: (M+)- ? ' 2 = for all j. Considering the limit Xq — ► 6 one similarly derives 
the condition (M+) J ' 4 = 0. In view of ({IB} one obtains more, namely 



whenever j = 2, 4 or A; = 2, 4. 



(54) 



Let us denote by M+' red the reduced 2x2 matrix obtained by omitting the vanishing 
rows and columns, i.e., 



Mt 



,red 



/(M+)M (M+) 1 - 3 
1 V (M+) 3 - 1 (M+) 3 -V 
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The condition $ 2 1 = for xq — * a means that the coefficient standing at (r 0a /2Y a exp(— 
i(a — l)0 0a ) vanishes. Using ((Mil we get 



sin(7r a) r(a) t 



2tt 2 



?1 sin(7ro;) r(a) 
7T 1 — a 



sin(7TQ;) r(a 



+ (M+)- 3 ^ "'^ (v^I) W §a- ll/3 -i(/3, a; *) ) = 0. 

7r 1 — a I 



This is equivalent to the couple of equations 



1 
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- (M+) 1 ' 1 (V=i) 2(1 " a) T^ 1|a _ a (a, /?; z) 



+ (M+) 1 ' 3 ( v^=I) M § Q -i^-i 09, «; *) 



- (M+) 3 ' 1 (v^) 2(1 ^ } T Q _i ja _i(a, /?; *) 



+ (M+) 3 ' 3 (V=i) 2 01 "Sa-i^-i^aj^) = 0. 

Analogously, another two equations are obtained when considering the limit Xo — » b, 
namely 



+ (M+) 3 ' 1 ( v ^) 2 ~ Q ~%_ w (a,/^) = 0, 
- (M+) 1 ' 3 ( v ^) 2(1 ^ ) 7fi-xfi-x{P, a- z) 

+ (M+) 1 ' 1 (^) 2 " M Vi, a -i(«,A2) = 0. 
The four equations can be jointly rewritten in the matrix form, 

M +,red 



(55) 



\_ ( W~zf 2Q T a _ liQ ,_ 1 (a,/3;z) -(v^)^ " P §/3-i, a -i(a,/?;^) 



2tt I -( v /37) 2 ° S^^^/?, a; z) (v^)^ ZP< J>-i,/3-i(/3, a; z) 



2-2/3 r 



It is straightforward to verify that the derived matrix M+ actually obeys conditions 
(jUJ) and (JIIJ). The former one follows from the equalities 



T MiV (q!, /3; z) = 7^ v (a, P; z), S M)1/ (a, /3; z) = 8^(0, /3; 2), 



and 



T^a, /5; z) = 7 Vtfi {a, p; z), S^(a, (3\ z) = S V)M (/3, a; z). 

The latter one follows from the form of P(z, w) given in (l53l) . In fact, (HE) and ifKB]) 
jointly imply 

(z-w)P(z,wY ed = (M+- red )" 1 - (M+' red ) _1 . 
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The other component of the Pauli operator, H , can be treated similarly. The 
boundary conditions read Q^ 1 = $2 = ( see (JH30- The condition 

= for 

means that the coefficient standing at (r 0a /2)~ 1+a exp( — i(a — l)9o a ) vanishes. Hence 

/•'(•'') <>• 

5 

or equivalently, (M") 5 '* 1 = 0. Similarly for a; — ► we derive that (M~) J ' 3 = 0, hence 
(M-) j > k = whenever j = 1, 3 or Jfe = 1, 3. (56) 

Set 

M -, red _ AM") 2 ' 2 (M") 2 > 4 \ 

The condition $2 = for xq — > a means that the coefficient standing at (r Qa /2) a exp(— ia9 
vanishes. Using (|56|) we get 



sin(7r a) T(l — a) 



2tt 2 



a 



/, a (*) + £/i 



;r 



,2a 



(M 



+ y',2 



sin(7r a) IY1 — a) 



7r a 

x (v /Z ^) ZQ 'J' Qia (a,/3; z) 
. + .,- 4 sin(7ra)r(l-a), , ^ a+l 3 , 



7r a 

This is equivalent to the couple of equations 



0. 



1 
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(M-) 2 ' 2 (v^i) 2Q T Q , a («,/5;z) + (M-) 2 ' 4 (v^) a+/3 S ai/3 (/3,a;z) 

For xq — > b one derives other two equations, 

{M;)^{V^z) 2P 7Ma ]Z ) + (M") 4 ' 2 (v 



1 

2~ix 



a+/3 



>,/?;*) 



-(M;) 2 - 4 (v^) 2 ^T^(/3,a^) + (M-) 2 - 2 (v^) Q ^§ fta (a,/5;z) 



a+/3 



0. 



Jointly the four equations mean that 



Ml 



.red 



1 

2^ 



(x/^) 2a T Q , a (a:,/3; z) -( v /3 ^) Q+P S/3, a (a,/3; z) 



(\/ = f) a+ a; z) (y^z)^ p 7p^{(3, a; 2;) 



2/3 r 



(57) 



Let us note that the inverted matrices on the RHS of ijHHI) and (IH7j) are actually 
well defined. This is because the matrices in question depend on z analytically in 
the domain C \ R + and tend exponentially fast to invertible diagonal matrices for 
Re y/—z — > +00 as one can easily deduce from the discussion of the formula (|32|) 
related to the convergence of the series (|31a|) and from the form of matrix entries (|38|) 
and (IH. 
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VII. Concluding remarks 



Having a formula for the Green function Sf(x, x ) it would be, of course, desirable 
to use it for a more detailed analysis of the Pauli operator, first of all for its spectral 
analysis. This aim would assume, however, a more detailed analysis of the functions 
&w,v((x, P', z) and 7 UjV (a, /3; z). In particular, it would be important to know what 
happens in the limit Re^/— ~z — > 0, i.e., when z approaches A G R+ from the upper 
or lower half-plane. Recall that both S UtV (a, /3; z) and T WjI/ (a, (3; z) are expressed as 
infinite series whose convergence is guaranteed for Re \f^z > 0. Our first attempts in 
this direction suggest that such an analysis might be rather complex and should be 
considered as an independent problem in its own right. 



Acknowledgements 

P. S. wishes to acknowledge gratefully a partial support from Grant No 201/01/0130 
of the Grant Agency of Czech Republic. V. A. G. is grateful to the Department of 
Mathematics of Czech Technical University for the warm hospitality and to INTAS 
(Grant No 00-257) and RFBR (Grant No 02-01-00804) for a financial support. The 
authors are indebted to the referee for helpful comments. 



28 



References 



1 Y. Aharonov, A. Bohm, "Significance of electromagnetic potentials in the quantum 
theory," Phys. Rev. 115, 485-491 (1959). 

2 Y. Aharonov, A. Casher, "Ground state of spin 1/2 charged particle in a two 
dimensional magnetic field," Phys. Rev. A 19, 2461-2462 (1979). 

3 L. Dabrowski, P. St'ovicek, "Aharonov-Bohm effect with 5-type interaction," 
J. Math. Phys. 39, 47-62 (1998). 

4 R. Adami, A. Teta, "On the Aharonov-Bohm Hamiltonian," Lett. Math. Phys. 43, 
43-53 (1998). 

5 P. St'ovicek, "The Green function for the two-solenoid Aharonov-Bohm effect," 
Phys. Lett. A 142, 5-10 (1989). 

6 S. N. M. Ruijsenaars, "The Aharonov-Bohm effect and scattering theory," Ann. 
Phys. (Leipzig) 146, 1-34 (1983). 

7 S. Albeverio, F. Gesztesy, R. H0egh-Krohn, H. Holden, "Point interactions in 
two dimensions. Basic properties, approximations and applications to solid state 
physics," J. reine angew. Math. 380, 87-107 (1987). 

8 P. Exner, P. Vytfas, P. St'ovicek, "Generalized boundary conditions for the 
Aharonov-Bohm effect combined with a homogeneous magnetic field," J. Math. 
Phys. 43, 2151-2168 (2002). 

9 O. Ogurisu, "Generalized boundary conditions of a spin-1/2 particle for the 
Aharonov-Bohm effect combined with a homogeneous magnetic field," mp_arc 
02-495. 

A. Arai, "Properties of the Dirac-Weyl operator with a strongly singular gauge 
potential," J. Math. Phys. 34, 915-935 (1993). 

1 V. A. Geyler, E. N. Grishanov, "Zero modes in a periodic system of Aharonov- 
Bohm solenoids," JETP Letters 75, 354-356 (2002). 

2 A. F. Nikiforov, V. B. Uvarov, Special Functions of Mathematical Physics 
(Birkhiiuser, Basel, 1988). 

3 J. Weidmann, Linear Operators in Hilbert Spaces (Springer- Verlag, New York, 
1980). 

4 N. I. Akhiezer, I. M. Glazman, Theory of Linear Operators in Hilbert Spaces (Pit- 
man, London, 1981). 

5 L. Dabrowski, H. Grosse, "On nonlocal point interactions in one, two and three 
dimensions," J. Math. Phys. 26, 2777-2780 (1985). 



29 



Figure captions 



FIGURE 1. Geometrical arrangement. Choice of the cuts L a , L b and choice of the 
angle variables 9 a , 9 b . 



FIGURE 2. Function i]) a ,a-i,\ from the deficiency subspace for the values of parame- 
ters a = 1/3, = 2/3, p = 1. 



FIGURE 3. Function ipb,/3,i from the deficiency subspace for the values of parameters 
a = 1/3, = 2/3, p = l. ' 
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Figure 1: Geometrical arrangement. Choice of the cuts L a , L b and choice of the angle 
variables 9 a , 9 b . 



31 



function ^a,«-i,/ 
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function if/b,p,i 
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